Graphene antidot lattice waveguides 



(N 

o 

(N 

> 
O 

OV 



o3 



O 

o 



> 
m 



(N 



X 



Jespcr Goor Pedersen 1 , Tue Gunst 2,3 , Troels Markussen 4 , and Thomas Garm Pedersen 1 ' 3 
1 Department of Physics and Nanotechnology, Aalborg University, Skjernvej 4A DK-9220 Aalborg East, Denmark 
2 Department of Micro- and Nanotechnology (DTU Nanotech), 
Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark 
3 Center for Nanostructured Graphene (CNG), Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark 
4 Center for Atomic-scale Materials Design (CAMD), Department of Physics, 



Technical University of Denmark, DK-2800 Kgs. 

(Dated: November 20, 2012) 



Lyngby, Denmark 



We introduce graphene antidot lattice waveguides: nanostructured graphene where a region of 
pristine graphene is sandwiched between regions of graphene antidot lattices. The band gap in the 
surrounding antidot lattices enable localized states to emerge in the central waveguide region. We 
model the waveguides via a position-dependent mass term in the Dirac approximation of graphene, 
and arrive at analytical results for the dispersion relation and spinor eigenstates of the localized 
waveguide modes. To include atomistic details we also use a tight-binding model, which is in 
excellent agreement with the analytical results. The waveguides resemble graphene nanoribbons, 
but without the particular properties of ribbons that emerge due to the details of the edge. We 
show that electrons can be guided through kinks without additional resistance and that transport 
through the waveguides is robust against structural disorder. 



PACS numbers: 73.22.-f, 72.80.Vp, 73.21.Hb, 73.21. Cd 



I. INTRODUCTION 



Graphene, the two-dimensional allotrope of carbon 
first isolated in 2004^*2 has emerged as a highly attrac- 
tive material for future electronic devices. Graphene 
has exceptional electronic properties, as demonstrated 
in its extremely high carrier mobility^ which even at 
room temperature is limited predominantly by impu- 
rity scattering^ Already, extremely fast graphene-based 
transistors have been realized,— and fabrication meth- 
ods have emerged, which allow for large-scale production 
of single-layered graphene^ One key element of future 
graphene-based electronics is the ability to localize car- 
riers in graphene wires, in order to facilitate electronic 
graphene networks. The most immediate way of realizing 
such wire structures is by cutting graphene into so-called 
graphene nanoribbons (GNRs)i 7 - Quantum confinement 
will in general induce a band gap in such structures, the 
magnitude of which scales with the inverse of the width 
of the GNR. However, the exact atomistic configuration 
of the edge of the ribbon greatly influences the mag- 
nitude of this gap, with particular edge configurations 
resulting in vanishing gaps or localized edge states^ 
An alternative to this is to rely on controlled genera- 
tion of extended defects in graphene to act as metallic 
wires^ Relying on results from two-dimensional elec- 
tron gases (2DEGs) formed at the interface of semicon- 
ductor heterostructures, one could also imagine applying 
electrostatic gating to define wire geometries. Graphene 
presents an interesting challenge in this regard via the 
phenomenon of Klein tunneling^ which makes it diffi- 
cult to achieve carrier localization in graphene via ordi- 
nary gating. Electrons impinging on a potential barrier 
at close to normal incidence are transmitted with van- 
ishing reflection, regardless of the height of the potential 



barrier. In spite of this, guiding of electrons via electro- 
static gating has been demonstrated experimentally^ al- 
beit with the caveat that guiding is restricted to a specific 
range of wave vectors for which Klein tunneling is negli- 
gible. In contrast, a so-called mass term in graphene pro- 
vides confinement that is a close analogue of gate-defined 
localization in an ordinary 2DEG. This term originates 
from a Dirac description of graphene, which emerges as a 
low-energy approximation of a tight-binding (TB) model 
on the honeycomb latticed Adding a diagonal term ±A 
to the resulting Hamiltonian, with the sign alternating 
between the two sublattices of graphene, the spectrum 
becomes that of gapped graphene, a semi-conductor with 
an energy gap of twice the mass term. This staggered 
potential is commonly denoted a mass term due to the 
analogy of the low-energy carriers in graphene with mass- 
less Dirac fermions, which acquire a mass via such a term 
in the Hamiltonian. 

In this paper, we propose realizing a waveguiding 
structure in graphene via graphene antidot lattices 
(GALs). GALs are nanostructured graphenes, which in 
their simplest description take the form of periodically 
perforated graphene structures.— GALs have recently 
been produced experimentally with both electron beam 
lithograph y 15 : 16 and block copolymer lithography] 17 i 18 
The periodic perforation induces a band gap in graphene, 
rendering it semi-conducting, and the resulting band 
structure closely resembles that of gapped graphene in 
the low-energy limit. GALs thus allow for the realiza- 
tion of a position-dependent mass term, as illustrated 
in Fig. [T] The idea is to sandwich a region of pristine 
graphene between two GAL regions, the band gaps of 
which define an energy range for which localized guided 
modes are expected to emerge in the central region. 
This idea is closely analogous to how photonic waveguid- 
ing is realized in photonic crystal structures^ We note 
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FIG. 1: (upper panel) Conceptual illustration of a graphene 
antidot lattice (GAL) waveguide. A central region of pristine 
graphene is surrounded by GAL regions, the band gaps of 
which confine states to the waveguide region. Translational 
symmetry is assumed in the y-direction, which is along the 
longitudinal direction of the waveguide, (lower panel) The 
geometry of a {7, 3}l z GAL waveguide. Black dots show the 
location of carbon atoms. Bloch boundary conditions are im- 
posed on all boundaries. The dashed lines illustrate the en- 
larged GAL unit cell, the width of which we denote by w. The 
lowest-energy waveguide mode at the F point, calculated via 
the tight-binding model, is illustrated with circles, the size of 
which shows the absolute value |^/>(a;,2/)| of the (real- valued) 
eigenstate, while the color indicates the sign. 



lattice basis vectors lie parallel to carbon-carbon bonds, 
to ensure that a band gap always exists for the GAL^ 
The lower panel of Fig. [1] illustrates the geometry of a 
{7, 3}| z waveguide. For simplicity we restrict the width 
of the waveguide to be an integer multiple of the width 
of the GAL unit cell. This is merely for computational 
convenience, and as we will demonstrate in this paper, 
simple scaling laws exist to predict the properties of more 
general widths of the waveguide. 

In this paper, two different methods will be employed 
to determine the waveguiding properties of GAL waveg- 
uides. We first consider a model based on the Dirac ap- 
proximation, including the influence of the GAL barri- 
ers via a position-dependent mass term. We will show 
that this model admits analytical solutions in certain 
limits, which is highly beneficial for determining the gen- 
eral dependence of the waveguide properties on the var- 
ious structural parameters. Furthermore, these results 
demonstrate clearly the unique properties of graphene 
waveguiding compared to quantum well structures de- 
fined in ordinary 2DEGs. To include the atomistic details 
of the structures we also consider a tight-binding model, 
which we use to calculate the transmission properties of 
the waveguides, taking into account potential structural 
disorder. 



II. DIRAC MODEL 



that other methods beside GALs have been proposed for 
achieving gaps in graphene via superlattices, such as, e.g., 
patterned hydrogenation^i or superlattices of boron ni- 
tride islands embedded in graphene^i Another alterna- 
tive is to sandwich graphene between hexagonal boron 
nitride, which is predicted to induce significant band gaps 
in graphene^ 

We denote a waveguide geometry as {L, H\™ ac , where 
L is the sidelength of the hexagonal unit cell of the sur- 
rounding GAL, while R is the radius of the perforations, 
both in units of the graphene lattice constant, a = 2.46 A. 
The width W of the waveguide is defined via TV = W/w, 
where w is the width of the enlarged GAL unit cell, as 
illustrated in the lower panel of Fig. [TJ The width of 
the waveguide is of course somewhat more ambiguous 
than in the case of GNRs, as we have no sharp edge 
defining the precise boundary between the region of pris- 
tine graphene and the surrounding antidot lattice. For 
denoting the geometries we simply take W as the dis- 
tance between the nearest edges of the two bordering 
antidot lattice unit cells. However, as we discuss be- 
low, the effective waveguide width is slightly larger than 
this. The width w of the enlarged GAL unit cell de- 
pends on the orientation of the waveguide with respect 
to the graphene lattice, which we indicate with the su- 
perscript, with 'zz' ('ac') denoting a waveguide with the 
longitudinal direction along the zigzag (armchair) orien- 
tation of the graphene lattice. Note that in both cases we 
choose the orientation of the GAL such that the supcr- 



A. Analytical derivation 

We first consider a simple model of a GAL waveg- 
uide based on a Dirac model of graphene. For graphene 
nanoribbons, the exact edge geometry can be included 
in a Dirac model via the boundary conditions of the 
spinor components in each valley^^ However, in our 
case, no atomically defined boundary exists between the 
central region and the bordering GAL regions and we 
thus adopt a model wherein the band gap, E g , of the con- 
fining GALs is included via a position-dependent mass 
term. We introduce dimensionless coordinates (x, 7) = 
2/W x (x,y) and dimensionless energies e = E/Eq, with 
Eq = 2%Vp/W ~ 1.278 eV-nmxVF -1 , assuming a Fermi 
velocity vf = 10 6 m/s of graphene. In these units, the 
eigenvalue problem for the spinor eigenstates reads 



m(x) ~i (d x - *<9 7 ) 
-i(d x + id 7 ) -m(x) 



<Mx>t) =e <Mx,t) 

<Mx,7)J |_<Mx,7) 

(1) 

where the dimensionless mass term m(x) = S = 
E g /(2Eo) for |x| > 1 and zero otherwise. This equation 
has been derived from the Dirac Hamiltonian, Hk, which 
emerges as a linearization of a TB model of graphene near 
the K point. We will discuss the differences between the 
two inequivalent K and K 1 valleys of graphene below. 
Note that the sign of the mass term is arbitrary, and 
that changing it has no physical significance, but merely 
results in an interchange of the two spinor components. 
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We stress, however, that the sign should be the same 
on both sides of the central, waveguiding region. The 
Hamiltonian commutes with the y-component of the mo- 
mentum operator and we thus take spinor components of 
the form <fi A (x,l) = /(x)e 4K7 and (j) B {x,l) = ff(x)e lK7 , 
where K = kW/2 is the dimensionless Bloch wave vec- 
tor along the longitudinal direction of the waveguide. 
As the mass term is piecewise constant, the equations 
for the spinor components can be decoupled in each re- 
gion of the waveguide structure. We look for bound 
states, and thus take e 2 < S 2 . The normalizable solu- 
tions for the first spinor component thus read /(%) = 
A±e ±l3x for |x| > 1 and f{\) = B cos(ax) + C sin(a%) 
for \x\ < 1 while the second component is given via 
9{x) = i{ K ~ 9 x )f(x)/(e + m{x))- Here, we have defined 



\/e 2 — k 2 and /3 



Vs 2 



The requirement 



of continuity of both spinor components at the bound- 
aries of the central waveguide region leads to a transcen- 
dental equation for the energies, V S 2 — a' 2 = a tan 2a, 
regardless of which valley is considered (see below). The 
solution of the problem is thus closely reminiscent of that 
of an ordinary one-dimensional square well potential, al- 
beit with the crucial difference that for graphenc E oc a 
rather than E oc a 2 , as a consequence of the linear disper- 
sion relation of graphene. Also, we note that, contrary 
to the case of the Schrodinger equation, the derivative 
of the eigenstate spinors need not be continuous at the 
boundary. 



B. Infinite mass limit 

The standard textbook graphical solution suggests 
that there are N bound states in the waveguide if 
(N — l)7r/2 < <5 < Ntt/2. In general, the energies and 
spinor components of these states will need to be deter- 
mined by numerical solution of the transcendental equa- 
tion. However, in the limit of an infinite mass term, 
S — > oo, the problem admits an analytical solution for 
the energies, 



C s ( K )= n / K 2 + ^(n+i) 2 , n = 0,1,2... (2) 



or, reverting to ordinary units, 



E 2 b (k) 



w 2 



(3) 



where Ef,(k) = hvpk is the bulk graphenc dispersion rela- 
tion and s = ±1. The dispersion relations of the waveg- 
uide modes thus resemble those of gapped grapheno 2 ^ 
with a mass term of A e s = (hvF^/W)(n + |). Se- 
ries expansion of the transcendental equation for the 
energies in the case of a finite mass term reveals that 
a ~ 71-/2(71 + 1/2)[1 + 1/(2S)]~ 1 so the results in the in- 
finite mass limit are expected to be valid when E g 3> 
2hvp/W ~ 1.278 eV-nmxW -1 . Including the leading 



order correction, the energies are given as 

m 2ft 2 « 2 7r 2 (n+i) 2 
E^{k) = EZ{k)-s- 



W 2 Jk 2 W 2 + 7r 2 (n+i) 2 



E a 



(4) 

The eigenstate spinors in the infinite mass limit are 
most easily determined by using the boundary condi- 
tions derived by Berry— which set a phase relationship 
between the spinor components at the edge of the waveg- 
uide. The spinor components for k — and e > can 
then be derived as 

f / \ / cos (n/2[n 4- l]x) for n zero or even, , . 
TAX) - \ S in(7r/2[n+ for n odd, 1 j 

i sin (7r/2[n + \]x) f° r 71 zero or even 



while 

9n(x) 



—i cos (7r/2[n + |]x) for n odd, 



for \x\ < 1. Here, we have omitted normalization con- 
stants. We note that because of the particle-hole sym- 
metry of graphene, the eigenstatcs should be normalized 
separately on each sublattice<21 The spinors for the e < 
eigenstatcs are given by exchanging /„ and g n . Inter- 
estingly, and in stark contrast to an ordinary infinite 
square well potential, this shows that in the limit of infi- 
nite mass, the charge density is evenly distributed within 
the waveguide, for states with vanishing wave vectors. 
We find that this holds true also in the more general case 
of a finite mass term, the main difference being that in 
this case the spinors extend slightly into the mass regions. 
Including a nonzero wave vector, the charge density is lo- 
calized predominantly on either the edges or the middle 
of the waveguide. 

The derivation above takes as its starting point a lin- 
earization of graphene near the K valley. The Dirac 
Hamiltonian near the inequivalent K' valley is such that 
the energy spectrum of the waveguide is the same for 
both valleys, while the eigenstate spinors are related by 
an interchange of the spinor components, i.e. <\)\ = (fig 
and (fin = ■ W e also note that to obtain the full wave 
function, a Bloch phase factor of e lK r or e lK ' r should be 
added to the spinors, with K and K' depending on the 
orientation of the graphene lattice. We will return to this 
point below, when we compare the results of the Dirac 
approximation with those obtained with a tight-binding 
model. 



III. TIGHT-BINDING MODEL 

To include the atomistic details of the waveguide struc- 
tures, and to clarify the validity of the results derived in 
the Dirac approximation, we employ a nearest-neighbor 
tight-binding approximation. This will also allow us to 
asses the significance of the orientation of the waveguide 
with respect to the graphene lattice. The TB model is 
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FIG. 2: (left) Band structure of the {7, 3}f GAL waveg- 
uide. The band structures are shown for the tight-binding 
model as well as the analytical infinite mass-limit results, 
E^ and numerical solution of the transcendental equation of 
the Dirac approximation, E ns . The shaded gray area shows 
the projected bands of the surrounding GAL regions. For 
comparison, the bulk graphene band structure, Et, is also 
shown. Note that A denotes the lattice constant of the waveg- 
uide, (right) Corresponding density of states for the TB 
model. Note the van Hove singularities characteristic of one- 
dimensional structures. 

parametrized via a hopping term t = — 3 cV between ir— 
orbitals, the on-site energy of which we set to zero. We 
ignore non-orthogonality of the 7r-orbitals. Fig. Q] illus- 
trates the geometry used for the TB model in the case 
of a {7, 3}| z GAL waveguide. We use periodic boundary 
conditions along the cc-axis as well as the y-axis. We 
have ensured that the results are converged with respect 
to the number of GAL unit cells included around the 
waveguide. The band gap of the GAL is quite well devel- 
oped even with just a few rows of antidots^ so including 
three GAL unit cells on each side of the waveguide usu- 
ally yields converged results. 

IV. RESULTS 

A. Dispersion relations 

In Fig. [51 we show the band structure of a {7, 3}| z 
GAL waveguide calculated using the TB model as well 
as the Dirac approximation. Only electron (E > 0) 
bands are shown. Both the TB and the Dirac model 
exhibit perfect electron- hole symmetry, so hole (E < 0) 
bands simply follow from a sign change. For the Dirac 
results we take the effective width of the waveguide to be 
W e s = (N+ \ )w, slightly larger than the definition used 
for denoting the waveguide geometries. Note that the 
wavevector is shown relative to the lattice constant of the 
GAL waveguide, which is A = 3La for the zigzag orienta- 



tion and A = \/3La for the armchair orientation. In the 
Figure, the shaded, gray region illustrates the projected 
bands of the GAL, which define the region below which 
localized waveguide states are expected to appear. This 
particular waveguide structure supports several localized 
states. Higher-lying band gaps also appear in the GAL, 
and we have confirmed that localized waveguide modes 
are also supported in these gaps. It is worth stressing 
that localized waveguide modes exist for all wave vec- 
tors in the first Brillouin zone. This is in contrast to the 
case of waveguides denned via electrostatic gating, where 
guided modes generally exist only for a limited range of 
wave vectors^ The dispersions of the waveguide states 
agree very well between the TB and the Dirac model, as 
long as the wave vector is not too near the Brillouin zone 
edges. The largest deviations between the two models oc- 
cur for energies close to the projected bands of the GAL, 
where coupling between the waveguide and the GAL 
states is pronounced. We note that the analytical result 
obtained in the infinite mass limit, E^ s is a very good 
approximation of the lowest waveguide mode. Includ- 
ing the first-order correction to E™ leads to near-perfect 
agreement with the full solution of the transcendental 
equation. Note that as illustrated in the derivation of 
the Dirac result, the waveguide dispersion relation quite 
closely resembles that of gapped graphene, which in turn 
is approximately the same as bulk graphene for wavevec- 
tors k 3> A c ft/ (hvp). This is evident in the figure, where 
for comparison we also include the bulk graphene dis- 
persion relation. The density of states (DOS) calculated 
using the TB model is shown in the right panel of Fig. [2] 
As expected from the Dirac approximation, the DOS 
closely resembles one-dimensional gapped graphene, i.e. 
g(E) (x Q(E - A cS )E/(hv Fy /E 2 - A^ ff ) with the van- 
Hove singularities characteristic of one-dimensionality 
clearly evident in the figure. 



B. Band gaps and effective masses 

In Fig. [3K, we show the energy of the lowest localized 
waveguide state at the T point as a function of the width 
of a {7,3}^ z ' ac ^ waveguide. Results are shown for the 
TB model as well as the analytical results obtained in 
the Dirac equation approach. Note that due to electron- 
hole symmetry, the band gap is twice this value. While 
the solution derived from the Dirac equation does not dis- 
tinguish between zigzag and armchair orientation of the 
waveguide, the TB model predicts that there are some 
differences between the two cases. To illustrate this, we 
show results obtained for waveguides oriented along the 
zigzag (ZZ) and armchair (AC) directions, respectively. 
The inset of the figure illustrates that both ZZ and AC 
orientations exhibit a clear 1/W dependence of the en- 
ergies, as predicted from the Dirac approximation, pro- 
vided the waveguide is sufficiently wide. We note that 
while differences do exist between AC and ZZ orienta- 
tions, these are rather small, and could be attributed 
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FIG. 3: (a) Energy of the lowest localized waveguide state at 
the r point, for the {7, 3} at family of GAL waveguides. The 
energy is shown as a function of the width of the waveguide. 
Results are shown for TB models (points) of waveguides ori- 
ented along the zigzag (ZZ) and armchair (AC) directions, 
respectively, as well as the analytical results obtained via the 
Dirac model in the infinite mass limit and including the first- 
order correction (black lines). See the legend in panel b. The 
inset illustrates the 1/W dependence of the energy for wide 
waveguides. For comparison, the thin colored lines show the 
energies for ZZ and AC GNRs, if certain edge dependencies 
are ignored (see text), (b) Corresponding effective masses 
in units of the free electron mass, m e . In both panels, note 
the close resemblance of the results obtained for AC and ZZ 
oriented waveguides. 



to a slightly different effective width of the waveguides in 
the two cases, due to the 7r/6 difference in the orientation 
of the surrounding GAL with respect to the waveguide. 
Indeed, results of the Dirac approximation fit both ori- 
entations quite well, especially when including the first- 
order correction. The results of the Dirac model can be 
made to fit even better if we take into account the fact 
that the effective width of the waveguide is likely to be 
somewhat larger than the definition we have used, see 
Fig. [I] Indeed, introducing the same effective width in 
the Dirac model as we did for Fig. [2] results in even better 
agreement with TB results. 

The absence of a well-defined edge means that the de- 
pendence of the properties of the guided modes on the 
width of the waveguide is much simpler than is the case 
for GNRs. In particular, GAL waveguides are always 
semiconducting, whereas in a nearest-neighbor TB model 
armchair GNRs alternate between metallic and semicon- 



ducting behavior depending on the exact width of the 
ribbon, while zigzag GNRs display dispersionless midgap 
states, localized on the edges£ In Fig. [3^,, we show the 
energies for ZZ and AC GNRs, calculated via the TB 
model. We stress that the ZZ GNRs also contain disper- 
sionless edge states at the Dirac point energy. Further- 
more, AC GNRs are metallic for widths W = (3j> — l)a, 
with p an integer £ To compare with the waveguide re- 
sults, we have included only semiconducting AC GNRs 
in the figure. With these modifications, there is quite 
good agreement between the energies of GNRs and the 
GAL waveguide structures. GAL waveguides thus resem- 
ble ribbons without the particulars resulting from edge 
effects. As such, we speculate that the electronic proper- 
ties of GAL waveguides may be easier to control experi- 
mentally. 

In Fig. we show the effective mass of the low- 
est waveguide state at the T point as a function of 
the width. From the analytical Dirac results in the 
limit of an infinite mass term, we find an effective mass 
m^j = hir(n + ^)/(vpW). Again, we note that there 
is excellent agreement between the results obtained in 
the Dirac approximation and those from a TB model, 
for both orientations of the waveguide. Including the 
first-order correction, results from the Dirac model are 
in near-perfect agreement with the TB model. While 
electrons in pristine graphene have vanishing effective 
masses, the appearance of an effective mass term in the 
waveguide structures results in non-zero, albeit still very 
small effective electron masses, which tend to zero in the 
limit of infinitely wide waveguides. 



C. Eigenstates 

To further compare the AC and ZZ waveguide orienta- 
tions, we show in Fig. @Ji the eigenstates corresponding 
to the lowest (positive) energy of {7, and {7, 

GAL waveguides, calculated at the T point. Note that 
these have approximately the same waveguide widths. 
The lower panels in the figure show the absolute square 
of the wave function, with the color indicating the sub- 
lattice. These results demonstrate a crucial difference 
between the AC and ZZ orientations, namely that while 
for the AC waveguide, the wave function is distributed 
evenly across the two sublattices, the ZZ waveguide ex- 
hibits pseudo-spin polarization, with the wave functions 
of the two sublattices localized predominantly on oppo- 
site edges of the waveguide. We note that the lowest 
energy is doubly degenerate, and that the second eigen- 
state (not shown) has the opposite pseudo-spin distri- 
bution. The upper panels in the figure show the inte- 
grated probability density, p(x) = J \ip(x,y)\ 2 dy along 
the transversal direction of the waveguide, with color in- 
dicating the sublattice. Note that despite the lack of a 
clearly defined edge, the probability densities very closely 
resembles those of GNRs^ The black line in the upper 
panels show the total probability density, if a broadening 
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FIG. 4: Wave functions of the localized waveguide modes corresponding to the (a) lowest and the (b) second-lowest (positive) 
energy at the F point of a (left) {7, and a (right) {7, 3}3 ac ' GAL waveguide. The lower panels in each case show the 

geometry, with carbon atoms indicated with black dots. Note that the actual computational cell includes additional GAL unit 
cells on each side of the central region. Superimposed on top of the geometry is the wave function, with the size of the circles 
indicating the absolute square of the 7r-orbital coefficient, while the color indicates sublattice. The upper panels show the 
integrated probability density, p(x) = f \tp(x, y)\ 2 dy. Red and blue circles indicate the densities on each sublattice. The black 
line shows the total density, when including a small broadening term. The dashed line shows the corresponding density at 
non-zero wavevectors, fcA = 7r/2. Note the rapid oscillations of the integrated density of the AC waveguide. 



of the order of the graphene lattice constant is included in 
order to account for the spatial extent of the 7r-orbitals. 
As predicted from the Dirac equation approach above, 
these results illustrate how the charge density is nearly 
uniformly distributed across the entire waveguide, also 
for the higher-lying states shown in Fig. HJd. This is in 
stark contrast to gate defined waveguides, which have 
wave functions more reminiscent of ordinary standing 
wave solutions^ In agreement with the Dirac results, 
we find that the uniform distribution only occurs at the 
r point. For non-zero wavevectors the density becomes 
localized predominantly at the edges of the waveguide for 
the lowest state, as illustrated with dashed lines in the 
figure for fcA = tt/2. In contrast to this, the densities of 
the second-lowest states tend to localize in the center of 
the waveguide as the wavevector is increased. 

To compare the wave functions with the spinor compo- 
nents derived via the Dirac equation, we first note that 
both the K and the K 1 points of graphene are folded 
onto the V point of the waveguide structure. We thus 
expect the eigenstates to resemble linear combinations 
of the eigenstates in both valleys. The differences be- 
tween the zigzag and armchair waveguides emerge due 
to the Bloch phase factors e lK ' r and e lK r , which differ 



depending on the orientation of the graphene lattice. For 
the zigzag orientation, the integrated probability density, 
p{x), becomes a simple linear combination of the eigen- 
states belonging to each valley, pa(x) oc \ Jk (x)+fK'{x)\ 2 
and pb{x) oc |<?jr(a:) + gK'{x)\ 2 . Using the expres- 
sions for the spinor components derived above, we find 
p(x) oc 1 ± sin(27r[n + ^]x/W), with the sign depending 
on the sublattice, which is in excellent agreement with 
the TB results. In contrast to this, because of the mix- 
ing of the valleys the probability densities of the arm- 
chair orientation exhibit a rapidly oscillating term, with 
a period 2tt/AK.— These rapid oscillations are clearly 
evident in the right panels of Fig. 2] As shown in the 
figure, these rapid oscillations are quickly washed out if 
a small amount of broadening is included. In this case, 
we recover the nearly uniform charge distribution within 
the waveguide predicted from the Dirac results. Finally, 
we note that the differences between armchair and zigzag 
oriented GAL waveguides are very similar to those seen 
between GNRs with armchair and zigzag edges.— In the 
case of GNRs, though, the difference emerges due to dif- 
ferent boundary conditions at the edge of the ribbon, 
which are not present in the case of GAL waveguides. 
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FIG. 5: Insets show the schematic transport setups for (a) a 
GAL waveguide and (b) a GNR, each connecting two semi- 
infinite graphene leads. Conductances for pristine (solid lines) 
and disordered systems (dashed lines) are shown for (a) GAL 
waveguides and (b) GNRs, in units of the conductance quan- 
tum Go- The straight black dashed lines show the graphene 
conductance. In the disordered systems, edge atoms have 
been randomly removed with a 5% probability. The length of 
the waveguides and GNRs are L = 89 nm. The shaded area 
in (a) indicates the energy range of the confined waveguide 
mode. 



D. Conductance 

Because the GAL waveguides have no clearly defined 
edge, one might wonder whether the guiding properties of 
the waveguides would be relatively robust to disorder. In- 
deed, the crucial ingredient is the existence of a band gap 
in the surrounding GAL. As this gap essentially occurs 
due to an averaging of the effect of the individual holes^ 
the emergence of a gap should be relatively robust to dis- 
order. A thorough investigation of disorder is beyond the 
scope of this paper, but as a preliminary study we model 
disorder by randomly removing atoms at the edges of the 
holes in the GAL. We consider a disordered GAL waveg- 
uide sandwiched between semi-infinite pristine graphene 
leads as illustrated in the inset of Fig. [5^- For com- 
parison we also consider the analogous system with the 



two graphene leads connected with a GNR having the 
same width as the corresponding waveguide (W — 4.5 nm 
for ZZ, W = 6.0 nm for AC), as illustrated in the inset 
of Fig. \5jp. We calculate the transmission through the 
waveguide and GNR using a recursive Green's function 
method^ with the lead-self energies determined using an 
iterative procedure.— The transmission is averaged over 
100 values of the transverse wave vector, and we fur- 
ther average over 10 samples with different realizations 
of the random disorder. To smear out the rapid oscilla- 
tions that occur due to interference between transmitted 
and reflected waves at the boundaries between the GAL 
waveguide and the graphene leads, we calculate the con- 
ductance at a finite temperature of T = 100 K. We con- 
sider a relatively high level of disorder, for which edge 
atoms are removed with a 5% probability. In Fig. [S£i, 
we show the conductance of disordered {7, 3}( zz ' ac ) GAL 
waveguides of length L = 89 nm (dashed lines) . For 
comparison, we also include the conductance of the pris- 
tine waveguides (solid lines). The shaded area indicates 
the energy range for the confined waveguide mode. In 
Fig. [5jj, we show the corresponding results for ZZ- and 
AC-GNRs. 

Focusing first on the conductances for the pristine sys- 
tems (solid lines), we observe that the ZZ- and AC ori- 
ented GAL waveguides have a similar conductance in the 
energy range of the waveguide mode. The high conduc- 
tances G ~ 1.5 Go show that there is a relatively good 
electronic coupling between the GAL waveguide and the 
graphene leads. The metallic ZZ-GNR is conducting at 
all energies, but the transmission saturates at a value 
close to unity. In the energy range of the waveguide 
mode, the GAL waveguides thus have a higher conduc- 
tance than both of the GNRs. Turning to the results for 
the disordered systems we observe that the GAL waveg- 
uides retain relatively high conductances. The ZZ and 
AC orientations show comparable reductions in the con- 
ductance due to disorder. The ZZ-GNR is less affected 
by disorder, while the AC-GNR conductance is signifi- 
cantly reduced except at the lowest energies. We con- 
clude that the GAL waveguides appear to be relatively 
robust against structural disorder and in general have 
higher conductances than the corresponding GNR sys- 
tems. 



E. Waveguide bends 

As mentioned earlier, the waveguides introduced in 
this paper arc closely analogous to photonic crystal 
waveguides. In such structures, light can be guided 
through bends in the waveguide with very little radi- 
ation loss.— Relying on this analogy, we expect GAL 
waveguides to show a similar robustness to kinks along 
the waveguide. To illustrate the localization of the elec- 
tronic state and the guiding properties of the GAL waveg- 
uide we show the local current through the waveguide in 
Fig. [SI Similarly to the results of Fig. [5J the waveguide 
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FIG. 6: Bond current through a 'kinked' {5, 2}J zz) GAL 
waveguide. The current is calculated at energy E — 0.25eV 
with a transmission of T = 1.9. The current is highly con- 
fined to the waveguide region and no additional reflections are 
observed due to the kinks. 



FIG. 7: Conductances of the 'kinked' waveguide shown in 
Fig. [6] (straight line) and of a straight waveguide of similar 
length (dashed line). The shaded area indicates the energy 
range of the confined waveguide mode of the straight waveg- 
uide. Note the nearly identical conductances of the two sys- 
tems. 



structure is connected on both sides to semi-infinite bulk 
graphene leads, not included in the figure. The left going 
bond current (per unit energy) in the presence of an in- 
finitesimal bias voltage is calculated from the left scatter- 
ing state spectral function, Al^j, and the hopping matrix 
elements Hij from the TB Hamiltonian. Between atom 
i and j the local current is Al^ Hij^ To visualize the 
current on the given scale the current running away from 
each atom was calculated and averaged over an applied 
mesh. The illustrated average current thus cannot be as- 
signed to the individual atoms anymore, which is the rea- 
son why current appears to occur within the holes in the 
figure. Since the /c-averaged transmission only changes 
slightly from the T-point result, we use the T-point scat- 
tering states. Fig. [S] clearly illustrates the confinement of 
the current to the waveguide region and the robustness 
against kinks. 

To further illustrate the strong guiding properties of 
GAL waveguides, we show in Fig. [7] the conductance 
of the waveguide bend illustrated in Fig. [51 A similar 
method was used as that for the results of Fig. [5] For 
comparison, we also show the corresponding conductance 
through a waveguide generated by omitting the kink in 
Fig. |5] and instead having the waveguide run straight 
through. These results show that while there are small 
differences between the two structures in the oscillations 
of the conductance, overall the introduction of a kink has 
almost no consequence on the conductance through the 
GAL waveguide. Very low reflection loss is thus intro- 
duced by the kink, despite the fact that the waveguide 
alternates between zigzag and armchair orientations. We 
note that we have found similar results for slightly dif- 
ferent waveguide structures. 



V. DISCUSSION AND SUMMARY 

The GAL waveguide systems studied in this paper rep- 
resent idealized structures, which may be challenging to 
realize experimentally due to the small hole sizes. How- 
ever, the applicability of the Dirac model allows for sim- 
ulations of arbitrarily large structures. Moreover, the 
Dirac model can equally well be applied to other gapped 
graphene systems, where the band gap is not induced 
through periodic holes, but e.g. via periodically absorbed 
hydrogen^ Although the Dirac and TB models applied 
in this work are very simple, previous studies of pure 
GALs have shown that both the Dirac- and TB models 
reproduce the trends obtained from more accurate den- 
sity functional theory calculations^ In any case, the con- 
cept of a GAL waveguide depends only on the existence 
of a band gap in the GAL region and not on the specific 
details, and we believe our simplified models capture the 
correct physics. 

In summary, we have introduced GAL waveguides. 
The band structures of GAL waveguides have been mod- 
eled with the Dirac model including a mass term, which 
is shown to be in excellent agreement with an atomistic 
tight-binding model. We have shown that GAL waveg- 
uides support modes which are highly confined to the 
waveguide region and are robust against structural dis- 
order and kinks in the waveguide. In transport calcu- 
lations, we find that GAL waveguides have higher con- 
ductances than corresponding graphene nanoribbons. A 
further advantage of the surrounding GAL may be that 
it will mechanically stabilize the structure and be able to 
carry some of the generated Joule heat away from the de- 
vice. GAL waveguides may thus be an attractable way of 
realizing electronic wires in integrated graphene circuits. 
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